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ABSTRACT 



The white dwarfs are promising laboratories for the study of cosmochronology 
and stellar evolution. Through observations of the pulsating white dwarfs, we can 
measure their internal structures and compositions, critical to understanding post main 
sequence evolution, along with their cooling rates, allowing us to calibrate their ages 
directly. The most important set of white dwarf variables to measure are the oldest 
of the pulsators, the cool DAVs, which have not previously been explored through 
asteroseismology due to their complexity and instability. Through a time-series 
photometry data set spanning ten years, we explore the pulsation spectrum of the 
cool DAV, G29-38 and find an underlying structure of 19 (not including multiplet 
components) normal-mode, probably £ = 1 pulsations amidst an abundance of time 
variability and linear combination modes. Modelling results are incomplete, but 
we suggest possible starting directions and discuss probable values for the stellar 
mass and hydrogen layer size. For the first time, we have made sense out of the 
complicated power spectra of a large-amplitude DA pulsator. We have shown its 
seemingly erratic set of observed frequencies can be understood in terms of a recurring 
set of normal-mode pulsations and their linear combinations. With this result, we 
have opened the interior secrets of the DAVs to future asteroseismological modelling, 
thereby joining the rest of the known white dwarf pulsators. 

Subject headings: stars:individual(G29-38) — stars: pulsations — stars: white dwarfs 

1. Introduction 

Due to their status as the end-products of most stellar evolution paths, the white dwarfs are 
among the oldest stars in the galaxy and therefore offer important clues about the universe around 
us. Studying their interiors will provide solid endpoints for stellar evolution, providing insights 



Department of Astronomy, University of Cape Town, Rondebosch 7700, Cape Province, South Africa 
Observatoire Midi-Pyrenees, 14 Avenue E. Belin, 31400 Toulouse, Prance 
Observatoire de Paris-Meudon, F-92195 Meudon, Principal Cedex, Prance 
'Department of Physics and Astronomy, University of Leicester, Leicester LEI 7RH, UK 
'PPARC Advanced Fellow 

Visiting Astronomer, Laboratorio Nacional de Astrofisica, CNPq, Brazil 
European Southern Observatory Casilla 19001 Santiago 19 Chile 
'JILA, University of Colorado, Box 440, Boulder, CO 80309 



- 4 - 



into nuclear reaction rates, mass loss mechanisms, and basic physical properties of matter in a 
variety of extreme conditions. Dating the oldest of these stars provides natural constraints on the 
age of our Galaxy, including possible independent measurements of the ages of the disk, the halo, 
and open and globular clusters. 

The key to exploiting the potential of the white dwarfs is buried beneath their relatively 
thin surface layers: we must discover their internal structure and composition to make them 
useful. Measuring the ages of the oldest white dwarfs, which have spent only a small fraction of 
their entire existence off the white dwarf cooling track, is as simple as measuring their cooling 
rates. Theoretical cooling rate models depend on the mass and structure of the white dwarf. 
Asteroseismology can both determine these physical parameters of a star and provide a means for 
directly measuring a star's cooling rate. 

The DAV (or ZZ Ceti) white dwarfs, with hydrogen-dominated spectra, are the coolest 
(f« 13, OOOK) and oldest of the four known classes of white dwarf nonradial g-mode pulsators.0 
As such, they are the most critical in answering the age question. Unfortunately, they have also 
proven the most difficult to understand. 

The techniques of asteroseismology work best when provided with an abundance of stable 



modes of oscillation. The first white dwarf successfully analyzed, the DOV, PG 1159-035, (Winget 



et al. 1991| ) provided over 100 such modes. The best analyzed DBV, GD 358 (|Winget et al. 19941) , 
obliged us similarly. Until now, however, the DAVs were not so willing to provide us with the 
necessary quantity of modes for easy analysis. 

The DAVs exhibit distinct trends with temperature: the hotter stars have lower-amplitude, 
shorter-period pulsations and the cooler ones have large-amplitude, longer-period pulsations. The 
hot DAVs have very few modes while the cooler DAVs have many, but most modes are unstable 
and seem to come and go without forming any obvious patterns of behavior. 

In his dissertation, Clemens (1993, 1994) found additional systematic properties in the hotter 
DAVs by looking at an ensemble of individual pulsators, thereby establishing a technique for 
applying asteroseismology on these objects which were thought to have too few modes for analysis. 
He was able to determine the individual properties of each of the hot DAVs he analyzed by seeing 
how it fitted in with group properties he discovered. In each star he found only a few modes, 
but when added together, they formed a much larger, coherent set of £=1 modes and a few £=2 
modes.i Clemens found a successive series of £=1 modes from k=l to k=6 even though no one 



Brown fc Gilliland 1994, for example, 



^The others three classes are the PNNVs, the DOVs, and the DBVs. See 
for a review of asteroseismology in a broader context. 

^The surface deformations on a pulsating white dwarf are usually modelled by the set of spherical harmonics, 
{Yl^{6,<l))} and an additional radial-node index, k. I is the total number of surface nodes and m is the number of 
surface nodes along a line of longitude. Most observed modes are £—1, with a smattering of £—2. There have been 
no reliable identifications of an £=3 (or higher) mode in any white dwarf star. 



- 5 - 



star had modes at each value of k. Each observed mode, however, was common to more than one 
star; when a star had an observed mode, it was always in one of the discerned groups. 

The implication of this result is that the overall structure of the hot DAVs must be very 
similar. Comparing the set of observed modes to theoretical models suggests the masses of the 
stars Clemens studied are quite near O.GMq and their hydrogen envelopes are near (where 

is the total mass of the star). The remarkable similarity of the hot DAVs supports the common 
asteroseismology credo that the pulsators are "otherwise normal stars," although the cool DAVs 
with their complex, variable power spectra remained, for the time-being, an enigma. The debate 
over "thick" versus "thin" Hydrogen layers, however, is far from over. See for example Fontaine 
et al. (1994), Fontaine and Wesemael (1997), and recent modelling results from Bradley (1997) for 
a recent summary. 

Here, we present results from an extensive study of one of the cool DAVs, G29-38. We find 
that although its power spectrum is not stable and changes quite dramatically from one observing 
season to another (with smaller changes within a season), we can still fit the observed set of 
pulsation modes with a set of predominantly i=l normal-mode pulsations. The key to this result 
was two-part: 1) obtain many seasons of data to observe a larger set of available modes, and 2) 
accurately identify and separate the linear-combinations from the more fundamental modes of 
oscillation. The linear-combinations are observed periodicities whose frequencies are sums and 
differences of those of other modes; we believe they arise mainly from non-linear effects in the 
system. 

There has been a great deal of previous interest in G29-38. Its variability was discovered 
by Shulov and Kopatskaya (1974) and confirmed by McGraw and Robinson (1975). Winget 
et al. (1990) observed it in 1988 and found a still unexplained phase variation of the dominant 
oscillation period, but did not explain the bulk of the star's pulsation properties. Zuckerman 
and Becklin (1987) reported a significant infrared excess in the star at wavelengths longer than 
2/im, the source of which is also still unknown, although orbiting dust is becoming the consensus 
(Zuckerman 1993, Koester, Provencal, & Shipman 1997). Barnbaum & Zuckerman (1992) report 



a possibly periodic radial velocity variation with an unknown source. Kleinman et al. (1994) used 
the phase timings of a stable pulsation mode at 284s to place severe limits on what kinds of 
orbital companions could be included in the system, showing the radial velocity variations cannot 
be orbital in origin. A recurring thread throughout most of these works is the variability and 
complexity of G29-38's power spectra. For the first time, we now propose a simple picture for 
G29-38's pulsations. 

Our results, combined with Clemens's work on the hotter DAVs, place all the DAVs 
alongside their other pulsating white dwarf cousins and open this crucial class of variables for 
asteroseismological analysis. We show there is a pattern behind the complex variable nature of 
these systems and, in the process, discover a different kind of data set is needed to solve these 
stars. A single set of observations, over a single observing season, no matter how well-resolved it 
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is, will not suffice. There must be well-sampled observations over many observing seasons. 



2. Techniques 

The white dwarfs are particularly rewarding objects for asteroseismology, the gains from 
which increase with the number of identified modes. (Modes are identified by specifying the values 
of the three integers, k, £, and m, along with z/ or P, the frequency or period.) Cepheids pulsate 



in one or two observed modes; the 5 Scuti star with the most known modes has around 20 (Breger 



1995). The white dwarfs, however, can pulsate in hundreds of observable modes. Each normal 
mode probes a slightly different region of the stellar interior, so having 100 modes is like having 
100 probes, all going to different depths and locations in the star's interior. 

The advantages of white dwarfs as asteroseismological laboratories quickly become their 
biggest disadvantage as well: they are very complicated. With so many modes active at the 
same time, we need extended data sets to be able to resolve and identify closely-spaced modes 
since the resolution of the Fourier transform (FT), which we use to search for periodicities in our 
reduced lightcurves, is proportional to the inverse of the time duration of the lightcurve. Even 
with a high-resolution FT, we can still have problems resolving closely-spaced modes if the alias 
peaks caused by observing gaps end up near real modes of the star. The FT finds an inevitable 
ambiguity in cycle count of each measured frequency due to these gaps. For nightly observations 
from a single-site, we therefore get alias peaks separated from the real ones by integral multiples 
of one cycle d~^. Unfortunately, this is near the typical white dwarf rotation rate which causes 
real modes to be present, separated by ^ ld~^ in frequency. 

In order to eliminate these troublesome ld~^ aliases, we strive to eliminate the one day gap 
in our data. To do this, we set up a network of collaborating astronomers around the globe, all 
observing the same star over the same time period with similar tools and observing techniques. 
This network, called the Whole Earth Telescope (WET: [Nather et al. 19"90| ), has been used 
quite successfully in the study of DO and DB variables. With the WET, we can obtain data 
uninterrupted by the daily rising of the sun and hence produce vastly improved FTs, with few 
aliases surrounding the real peaks. 

Armed with a nearly alias-free transform, the goal of asteroseismology is to match each 
observed frequency with a unique value of k, i, and m. If a full set of modes (say all the possible 
m values for at least two values of i over a consecutive series of k values) are present in the star, 
the job is relatively easy. If not, then clues must be taken whenever available and pieced together 
for a consistent final picture. The clues involve the spacings between modes of same k and i but 
different m, and those between modes of same i and m, but different k. 

As the number of radial nodes {k) increases, the frequency of a g-mode decreases (there is less 
of a restoring force since the wavelength of the oscillation decreases, meaning less mass is present 
to supply the gravitational restoring force). In the asymptotic (high k) limit, the modes with same 
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a are equally spaced in period. The periods P, of such modes are given by (see, for example Unno 
et al. 19891 ): 

A;An 

+ constant (1) 



^^w+^ 

where AH is a constant related to the period spacing and the additive constant is small. AH 
itself is primarily a function of the mass of the star and is truly constant only for stars of uniform 
composition. Adding compositionally stratified layers (white dwarf stars have a very high gravity 
which separates and stratifies the constituent elements), or any other radial discontinuity, to a 
model star makes the value of AH different for each mode, although the mean remains a good 
measure of the total stellar mass. The deviations from uniformity of AH measure the layering 
present in the star. This effect is called mode-trapping as modes with nodes near the composition 
transition boundaries will have their periods shifted slightly so the nodes correspond more closely 
to the transition discontinuities ( [Winget, Van Horn, fc Hansen 1981 , Brassard et al. 1992). 



As long as the spherical harmonics are valid representations of the observed surface distortions 
(they are as valid as our assumption of spherical symmetry), their underlying symmetry implies 
that the period of each mode depends only on A;, and the total number of surface nodal lines 
and not on how the nodes are arranged on the surface (i.e. , m). When the underlying symmetry 
is broken, however, the observed periods become a function of m as well. Rotation breaks the 
symmetry and splits each mode of a given k and £ into a multiplet of modes with 2^ + l components 
with m running from to i. 

For slow rotation, the frequency difference (Ao") for each m mode is given, to first order, by: 



Act = m(l - Ct^k)^ (2) 

where is the constant stellar rotation frequency and Ci^k is, in general, a complicated function 
of the star's density and modal displacements which in the asymptotic (large h) limit, approaches 
the value 1/^(^+1). Thus, modes of the same k and but different m will be uniformly spaced in 
frequency and, barring radially differential rotation, all modes with the same i will have the same 
frequency spacings. If radially differential rotation is present, which changes the effective value of 
Q for each mode, we will see a corresponding systematic pattern in the m-splittings as a function 
of k. Since the white dwarf rotation periods so far measured with the WET are near one day, 
we expect the i=\ m-spacing to be on the order of 6^Hz, much smaller than the typical period 
spacing, AH, which is close to 150//Hz in the region of main power. Since these rotational spacings 
are small, we expect to find closely-spaced triplets for 1=1 modes and quintuplets for i=2. 

Having accurately isolated and measured the frequencies of the pulsations, with a long 
timebase of observations we can actually observe stellar evolution in progress. The period of a 
mode changes slightly as the star cools and/or contracts. The contraction process (presumably 
dominating only in the hottest of the pre-white dwarf pulsators, since white dwarfs cool at 
essentially constant radius) decreases the period with time while cooling increases it. Since the 
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rates of change (or Ps) are very small (of order 10~^^,s/,s for an average DAV), this is a very 
difficult measurement, but an extremely important one. Once we can measure the rate at which 
a star is cooling (Winget et al. 1991, Kepler 1993[ ), we can empirically calibrate the white dwarf 
cooling curve, and hence the luminosity function ( [Liebert, Dahn fc Monet 1988 ) to measure 
directly the age of the white dwarfs (Winget et al. 1987). 



3. Observations 

The WET is the ideal tool for analysis of stars with stable pulsation spectra. For most DAVs, 
however, it is not enough. During any given season of observations, no matter how well resolved 
the power spectrum of any one star, only a few independent pulsation modes are present; the 
rest of the frequencies found in the power spectrum are linear combinations of existing modes. 
If, however, we observe the same star over many seasons, we can add a few new modes each year 
and slowly build them up to see a larger set of possible modes. This approach will work only if 
there is a stable underlying pattern to the star's pulsations: that is, if it picks and chooses a few 
modes each season from a pre-defined, limited set of possibilities. We found this to be the case 
with G29-38. 

Presented here are over 1100 hours of time-series photometry (we call the technique temporal 
spectroscopy) on G29-38, representing the results of three campaign observations (two WET runs 
and a double-site venture between SAAO and McDonald) and many years of intense single-site 
coverage. They span the 10 years from 1985 to 1994 with data every year except 1986 and 1987. 
Rather than present the complete lengthy table of observations here, which can be found in its 
entirety in Kleinman (1995) and nearly complete in Kleinman et al. (1994), we list in Table |l| only 
the sites and telescopes which have contributed data, and in Table |2| the additional observations 
used in this work that were not listed in the Kleinman et al. (1994) paper. 
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Table 1. Data-providing sites. 



Location 


Telescope 


CTIO 


1.5m 


Itajuba, LNA 


1.6m 


KPNO 


1.3m 


La Palma (INT) 


2.5m 


Maidanak 


1.0m 


Mauna Kea (Air Force) 


24" 


Mauna Kea (CFHT) 


3.6m 


McDonald 


30" 


McDonald 


36" 


McDonald 


82" 


OHP 


1.93m 


SAAO 


30" 


SAAO 


40" 


SAAO 


74" 


Siding Spring 


24" 


Siding Spring 


40" 
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Table 2. Additional observations used in this work that are not hsted in 

Kleinman et al. (1994). 



Telescope 


Run Name 


Date 
(UT) 




Start 
(UT) 


SAAO 74" 




198^ Alio- 


» 

O 




S A AO 74" 


qSfiflfi 
oOUUU 


1Q8^ Alio- 


in 






s361 


±C>00 iT. LLg 


1 3 

-L (J 




SAAO sn" 

kJiT.iT.vy ou 


s361 8 




1 4 


9n-93-4n 


SAAO SD" 


s3621 


1985 A IIP- 


1 5 

X (y 


20-94-00 


SAAO SD" 

kJiT.iT.V-/ *J\J 


s3624 


1985 Aup- 




1 9-51 -20 


SAAO W 

kJiT.iT.V-/ OVJ 


s3fi28 


1985 Ane- 


17 

J. 1 


20-14-20 


SAAO W 

kJiT.iT.\_/ OVJ 




1985 Aug 


19 


20-42-20 


SAAO 4f)" 

kJiT.iT.vy 


s3634 


1985 Aug 


20 


91 ■05-1 6 


SAAO 40" 

kJiT.iT.vy '±VJ 


s3638 


1985 Aug 


21 


'22■2^ -45 


MrDnnald W 


r3084 


1985 Aug 


22 


6-37-26 




r308B 


1985 Aug 


23 


7-43-21 


SAAO 40" 

kJiT.iT.Vy t:V7 


s3fi42 


1985 Aug 


23 


20-54-18 


SAAO 40" 

kJiT..iT.Vy AW 


s3645 


1985 Aug 


24 


20-59-43 


QA AO 4n" 

kJiT.iT.vy t:VJ 


s3fi47 


1985 Aug 


25 


20-06-20 

.\J\J .^\J 


MrDnnalrl Sfi" 


r3086 


1985 Aug 


26 


7-19-19 


SAAO 40" 

kJiT.iT.V_/ 'iU 




1985 Aug 


26 




SAAO SO" 

kJiT.iT.Vy OVJ 




1985 Sep 


10 


18:44:20 


SAAO 30" 

kJiT..iT.Vy fJ\J 


s3655 


1985 Sep 13 


17:42:00 


SAAO 30" 

kJiT..rT.vy 


s3fi56 


1985 Sep 14 


17:43:00 


SAAO 30" 

kJiT.iT.Vy (JW 


s3658 


1985 Sep 15 


0:42:40 


SAAO 30" 

kJiT.JT.Vy OVJ 


s3660 


1985 Sep 15 


17- 39 -20 


SAAO 30" 

kJiT.JT.Vy *J\J 


s3663 


1985 Sep 16 


17- 39 -20 


MrDonalH 82" 


r3088 


1985 Oct 22 


2-58-37 


MrDonalH 82" 


r3094 


1985 Oct 31 


3-02-00 


McDonald 82" 


r3095 


1985 Nov 1 


1:56:30 


McDonald 36" 


sjk-0264 


1993 Jul 21 


9:16:30 


McDonald 36" 


tkw-0034 


1993 Aug 


11 


8:23:00 


McDonald 36" 


tkw-0040 


1993 Aug 


16 


4:16:00 


McDonald 36" 


sjk-0265 


1993 Sep 14 


2:04:00 


McDonald 36" 


sjk-0266 


1993 Sep 14 


9:24:30 


McDonald 36" 


sjk-0267 


1993 Sep 15 


1:53:30 


McDonald 36" 


sjk-0268 


1993 Sep 16 


2:03:30 


McDonald 36" 


sjk-0269 


1993 Sep 17 


2:34:00 


McDonald 82" 


sjk-0270 


1993 Sep 18 


3:52:30 


McDonald 82" 


sjk-0276 


1993 Sep 19 


3:39:30 
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The data were all reduced as described in Nather et al. (1990), transforming the raw 
lightcurves to extinction corrected relative intensity measurements. The raw data consist of either 
two or three lightcurves, measuring the variable star, the comparison star, and in the case of 
three-channel instruments, the sky as well. With two-channel instruments, we occasionally move 
the telescope off the target and comparison stars to sample sky. The instruments used for most of 
these observations are described by Kleinman et al. (1996). 

4. Results 

After reducing each data set, typically a week or two in length, we calculated its FT and show 
them schematically in Figure ||. Each line represents a periodicity identified at that period. Since 
the dynamic range in the system is large, we have ignored amplitude information and plotted 
every periodicity with a line of the same height. Each panel is labelled for the month and year of 
the data, with the exception of X2N88 and X8S92, which are data from the second WET run in 
November, 1988, and the eighth WET run in September, 1992, respectively. The bottom row of 
the plot, labelled Sum, contains all the modes in the panels above, collapsed into one. 

While there is a wealth of information in this plot, the most striking feature is the near 
continuum of modes shown in the Sum panel. Were all these modes simply i=l and even i=2 
modes, the mode density would not be nearly so great and there would be distinct gaps between 
modes. 

The most obvious explanation for this result is the known abundance of linear combination 
frequencies which were not removed in the plot. Our ability to identify (and hence remove) 
the linear combinations depends greatly on the signal to noise ratio and the resolution of each 
transform. To help avoid uncertain and incorrect identifications, we now restrict the analysis to 
the best data sets available — one per year: Aug85, a two-site campaign by South Africa and 
McDonald; X2N88, the first WET campaign; Sep89, a 20-night data set from McDonald; Oct90 
and Sep93, slightly smaller single-site data sets from McDonald; and X8S92, the second G29-38 
WET campaign. 

Identifying the combination frequencies is a difficult task. First, we must identify the three 
frequencies that form the combination, then decide which are the "real modes" and which is the 
combination. That is, if we find three signals with frequencies. A, B, and C such that A-|-B=C, 
we need to decide if A is the difference frequency of two modes B and C, if B is the difference 
frequency of A and C, or if C is the sum frequency of A and B. To identify the combinations, we 
wrote a simple computer program that goes through a list of identified frequencies from the power 
spectra. In developing this code, we found it identified every combination we found by hand, plus 
a few we did not. In addition, while we often stopped after finding one possible combination 
that contained a given mode, the computer code found all possible combinations, often finding 
some more exact than those we had originally identified. The code takes into account possible 
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Table 2- 


-Continued 




Telescope 


Run Name 


Date 


Start 






(UT) 


(UT) 


McDonald 82" 


sjk-0277 


1993 Sep 20 


3:39:00 


McDonald 82" 


sjk-0278 


1993 Sep 20 


5:05:30 


TV /T 1 1 r)0?5 

McDonald 82 


sjk-0279 


1993 Sep 20 


10:32:30 


McDonald 82" 


sjk-0281 


1993 Sep 21 


4:38:30 


McDonald 36" 


sjk-0282 


1993 Nov 4 


1:11:30 


McDonald 36" 


sjk-0283 


1993 Nov 5 


1:12:30 


McDonald 36" 


sjk-0285 


1993 Nov 6 


0:54:00 


McDonald 36" 


sjk-0287 


1993 Nov 7 


0:55:30 


McDonald 36" 


sjk-0288 


1993 Nov 8 


2:17:00 


McDonald 36" 


sjk-0289 


1993 Nov 8 


5:45:00 


McDonald 36" 


sjk-0293 


1993 Dec 14 


1:57:30 


McDonald 36" 


sjk-0295 


1993 Dec 15 


0:38:00 


McDonald 36" 


sjk-0297 


1993 Dec 16 


2:05:30 


McDonald 36" 


sjk-0299 


1993 Dec 17 


1:01:30 


McDonald 36" 


sjk-0302 


1993 Dec 18 


0:49:00 


Siding Spring 24" 


sjk-0363 


1994 May 14 


18:44:00 


Siding Spring 24" 


sjk-0368 


1994 May 15 


18:44:30 


Siding Spring 24" 


sjk-0373 


1994 May 16 


18:43:30 
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Fig. 1.— 



Schematic diagram of G29-38's periodicities for the entire data set. 
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mis- identification of the frequencies if there are dominant ahases present, and apphes a selectable 
equality criteria in determining each match. 

Once we had identified the combinations, we decided which were real and which were 
combinations using the following guidelines: relative mode amplitudes (combinations generally 
have smaller amplitudes), number of combinations with each mode (combinations with combination 
modes are less likely than first order combinations), the existence of multiplet structure (two 
multiplets adding together produce distinctly different structure in the combination), and the way 
modes and their combinations disappear and recur in the various data sets (a parent mode can 
appear without its combinations, but not vice versa). Still, though, there were a couple occasions 
when we could not decide without resort to a model that predicts normal modes. In these cases, 
we assumed modes with periods near the expected groups were real and not combinations. We 
will gladly make available the list of frequencies or transforms to anyone who wishes to seriously 
explore this procedure in more detail. The FTs and mode lists are included in Kleinman (1995) as 
well. 

Figure ^ is the schematic period diagram for this data subset, after removing the identified 
linear combination frequencies. This new schematic period diagram is substantially cleaner 
and shows the mode groupings we would expect for normal-mode pulsations. The roughly 
equally-spaced groups seen in the Sum row suggest a mean period spacing of roughly 50s, 



consistent with models of O.GMq DAVs i=l spacings ( Bradley 1996 ). If these are all i=l modes, 
we see a nearly complete set of them from 110s to 900s with two more at longer periods. With few 
exceptions, the groups are tight and have distinct gaps between them. 

As Figure |l] readily shows, the power spectrum of G29-38 changes dramatically from year 
to year. This is even more obvious in Figure ^ which shows a portion of the FT, retaining real 
amplitudes, for several of our best data sets. It appears, however, to make its most dramatic 
changes when the star is behind the sun and we cannot watch it. (This certainly sets an upper 
limit on the timescale of change, something that might prove useful once we understand why it 
makes such drastic changes.) While this predisposition means we do not get to watch the star 
change, the advantage is that the FT of each individual year provides a relatively stable set of 
modes that may be used for asteroseismological analysis. 



5. Mode Discussion 

As discussed earlier, we expect g-modes of identical £ values to be equally spaced in period in 
the asymptotic (high-A;) limit. We cannot, however, now go and search for strictly uniform period 
spacing, because we are not in asymptopia here. We do not expect to find a strictly uniform 
period spacing since we have identified modes possibly starting at k=l, which is certainly not in 
the high-A: limit. Instead, we expect to see significant departures until we get to the higher-A: 
modes (our experiences shows A:=10 is a good rule of thumb). We also expect deviations on the 
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Schematic diagram of G29-38's periodicities minus the hnear combination modes. 




Fig. 3. — A portion of G29-38's FT for several of our best seasons. The large amplitude changes 
are obvious. 
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order of 10s or so due to mode-trapping effects. In some cases, however, models of Bradley (1996) 
show deviations of as much as 20-30s. 

We tried many different methods to quantitatively search for equal period spacings in this 
observed set of modes while simultaneously allowing for the inevitable departures from strictly 
uniform spacing, but did not succeed in finding any significant results. The fault may not, 
however, be in the data: the same tests failed to find significant period spacings in some of the 
model-derived test data as well. We therefore believe that the inability to statistically determine 
a significant mean period spacing does NOT affect (either pro or con) the model of single-^ 
pulsations and continue under the assumption that the approximate equal period spacing seen in 
these modes is significant and that most likely, these are a series of successive k, £=1 g-modes. 
While the goal of studying this star, and others like it, must be to uniquely determine the k, i, 
and m for each mode so we can proceed to explore the white dwarf interiors, were we to go no 
further in the analysis than where we are now, we would be concluding with the discovery of a 
a pattern to G29-38's seemingly haphazard behavior: a recurring set of stable modes. With its 
now understood wealth of normal mode oscillations, G29-38 has become the DAV of choice for 
asteroseismology . We could now commence additional observing programs to help determine 
i values including perhaps, the spectroscopic methods described by Robinson et al. (1995). In 
addition, we now know how to make use of the entire set of easy-to-observe, large-amplitude DA 
pulsators, to get a similar set of information. 

However, since we are compelled by the observed period spacing, we will continue the analysis 
under the assumption that each group is a different k, £=1 mode. In figure ^, we provide a running 
k assignment for each group. We will later discuss the uncertainties of this assignment, but for 
now, they serve as references to make the ensuing discussion easier. 

The data from September, 1993 have the most modes and nicely reproduce the ~ 50s spacing 
seen in the Sum row earlier. This season is unique in that there appear to be six consecutive 
overtones (A;s) in one period range. The other seasons also show roughly equally spaced groups, 
but are often missing one or more A:s in between each observed mode. The Sum row shows all the 
gaps (but one) are filled and we see what appears to be a set of same-^ modes from k=l to k=17. 
Almost half of the observed modes repeat at least once in the data set and four are present in four 
of the five data sets. We now have strong evidence for a series of successive- A;, same-i modes. 

These same-£ modes fit well with an i=l model and its expected period spacing with the 
possible exceptions of the two modes near k=17, although they could also be k=17 and A;=18 
without much problem. These two modes, the 894s mode in Sep93 and the 915s mode in X2N88, 
are separated by 21s, a little too close, but not completely impossible, to be different ks (17 and 
18), same-i pulled closer by mode trapping. Models of Brassard et al. (1992) and Bradley (1993) 
with hydrogen layer masses near 10~^''M^ do, for example, show the strong mode trapping this 
identification implies near these periods, although this match requires a rather limited parameter 
set. In frequency, the separation is 26 /xHz, too large to easily fit the frequency spacings of the 
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known multiplets. We cannot yet completely settle this ambiguity, but are relieved to know it 
does not affect the rest of the pattern. 

5.1. Identification Ambiguities 

As the two footnotes in our list of observed periods, Table ^ suggest, there are still some 
ambiguities in our mode identification assignments. The two most critical uncertainties, as luck 
would have it, are for the two modes most critical in model- matching: the lowest-^ modes. The 
k=l mode seen in the X2N88 data set has two combinations: plus and minus the k=10 mode. 
There are no other combinations with these modes, we don't see them in any other year, and their 
amplitudes are all quite similar, thus it is impossible with this data to determine which mode is 
not a combination mode. The most likely choice is that we are seeing the k=l mode at 110s and 
its sum and difference with the k=10 mode. However, we could also be seeing the k=l mode at 
134s with two sums with the k=10 mode. It is unlikely, given the preference for observable sums 
over differences, that both the 110s and 134s are differences with the higher frequency k=l mode 
at 93s. 

The k=2 mode identification seen in the Sep89 data set is also uncertain, but this time not 
due to the presence of combinations, but rather the lack of them. This mode is of extremely 
low-amplitude and is not well-resolved in the FT. Because of the very high-quality data in this 
data set, the power seen in the region is probably significant, but we cannot say for certain where 
the mode is or even if it is stable over the course of the run. Like the k=l mode, we don't see this 
mode in any other data set. Placing the mode at the highest peak in the region and identifying it 
as 1=2 fits in well with the expected mode pattern, but is highly uncertain. 

Also in the X2N88 data is a messy region of excess power near 2105/iHz, its largest-amplitude 
component. We also see a series of three combination sums of this region (and the 2105)uHz peak) 
with the large-amplitude A;=10 mode. One of either the 2105/xHz peak, or its sums with the k=W 
mode, therefore, must be a real oscillation mode, but we cannot decide which it is. No matter 
which mode it is, however, it does not appear to fit in well with the established pattern of i=l 
modes and their associated period spacings. 

In the Sep89 data, we see two modes at 1986 and 2020 /iHz, near what we have been calling 
the A;=8 region, but also significantly different from the power seen there in other years. As seen 
in the next section, this mode shows strange multiplet structure and the modes seen here could be 
some artifact of that. Or, they could be something else entirely; we cannot yet say. 

Also in the Sep89 data are two modes at 2747/iHz and its sum with the 1986^Hz mode. The 
2747/iHz mode is about 3 times larger in amplitude than the sum mode, and thus is likely the 
non-combination mode, but we have no other evidence to support this. Neither mode fits in well 
with our i=l spacings. 
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We suspect therefore, that we do have some non-£=l modes in our data, but the lack of 
repeatabihty and consistently low-amplitude makes it difficult to say with certainty. If we accept 
the i=l assignments for the k=l, k=2, and perhaps even the k=8 modes (1986 and 2020 /xHz) 
of Sep89, then we have only two modes, the 2105//Hz group of X2N88 and the 2747/iHz pair of 
Sep89 to explain some other way or with some other i. 

6. Multiplets 

While we have uncovered strong suggestions for the non-combination modes being all i=l, we 
would like additional proof. Ordinarily, rotational splitting of different m values for a given i and 
A; is a valuable aid in assigning ^-values. Unfortunately, we do not find many multiplets in this 
data set. Those that are present do not behave quite as theory predicts. Since a good picture is 
often worth a few pages of text, we have schematically plotted in Figure § the resolved multiplets. 
This time, the height of each arrow is proportional to the amplitude of the mode (although the 
scale changes for each multiplet series) and the shorter line segments are not real modes, but 
simply represent the average period of the two modes on either side of it and is meant to suggest 
the location of the presumed missing m=0 mode. 

Only twice do we see identical spacings within a given data set. In the X2N88 data, the 
500s and 615s multiplets both have a 8.5/iHz spacing. In the Sep93 data, the 400s and 612s 
multiplets both have an ~ 4.7/iHz spacing. The 400s multiplet always has the smaller spacing; the 
500s multiplet, the larger; and the 600s modes alternately fall into both camps. Also note the 400s 
multiplet spacing is not constant, but varies (perhaps periodically) with time. All these results 
are quite interesting, and likely to be clues to the complicated structures in this star, but seem 
to be largely irrelevant to the overall mode picture. That is not to say there is no information 
in these changes, however; they may very well be understood by some nonlinear effects such as 
those discussed by Buchler, Goupil, and Serre (1995). We have only recently begun to acquire the 
kind of data that could comfortably test such theories. The unfortunate thing here is they do not 
directly yield information about the £-value of our modes, but seem to be even more mysterious if 
we invoke multiple-fe than if we don't. Thus they support, but do not demand, our £=1 model. 

7. Analysis 

We have now discovered a fairly regular pattern of recurring modes in the power spectra of 
G29-38 and have suggested they may all form a pattern of successive k, i=l g-modes. Going 
just this far is a very important step to uncovering the asteroseismological secrets hidden inside 
these stars; we never before knew if the modes we were seeing were sustained, normal-modes of 
oscillation or something much more fleeting and chaotic. This is just the first step, however. What 
remains now is to explore models of DAVs, in attempts to tune the input physics to the observed 
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set of frequencies. Such modelling work is a large task and is best allotted a separate work for 
the detailed, careful analysis and discussion which must inevitably occur. However, here we can 
address the direction such an effort might take. We do not intend these results to be concrete, but 
they should help point the modelling investigation in a productive direction. 

The goals of the detailed investigation will, of course, be to constrain or determine the mass 
and internal structure of G29-38. The asteroseismologically-determined estimates of both the 
total stellar and hydrogen layer masses will be particularly important in helping to both calibrate 
the asteroseismological methods with other mass determinations and in checking models of white 
dwarf formation and evolution which set limits on the amount of hydrogen in the DAVs. We also 
want to see how G29-38 fits in with the other DAVs, both the hot and cool ones: does it fit the 
period structure of the hot DAVs found by Clemens (1993, 1994) and does the period spectra 
of G29-38 look like that of the other cool DAVs, to the extent we can tell? We must defer the 
full answers to these questions to papers now in progress (however, see [Kleinman 1995 , Dolez & 
Kleinman 1997| ), but will address some of the key issues in a preliminary way here. 



If this observed set of modes is a set of single-^ pulsations, then the mean period is an 
indicator of mass and the deviations from the mean describe the layered internal structure of the 
star. Using forward differences (thus defining AP), we have plotted in Figure ^ a AP vs. P, or 
period spacing, diagram for these data. With the possible exception of the two modes near k=17, 
there is no reason to suspect any of these modes are a different i (we later treat these two k=17 
modes as k=n and A:=18). Based on the evidence already discussed and previous asteroseismology 
from WET observations and the work of Clemens with the hot DAVs, it seems unlikely that 
the majority of the modes are anything other than i=l. However, since it is possible we have 
inadvertently included a mode or two of different i into our identifications, we have plotted some 
likely alternatives to the the i=l only model by dotted lines and circles. Where we have plotted 
such alternatives, the only i=l interpretation is plotted with thinner lines and slightly smaller 
filled circles. The dotted squares represent the average spacing for two modes with an assumed 
missing mode between them. The bottom panel of the plot is an enlarged plot of the sum row of 
Figure || with arrows showing which modes were used to calculate the differences shown above it. 
Above this box are 1=1 k assignments. Each point has an additional uncertainty of a few seconds 
since we cannot be certain of the m values. In the worst case, the shift can be as much as 14s for 
periods near 900s, and 3s for periods around 400s, assuming the largest splitting of 8.5/iHz. With 
the smaller 4.7^Hz splitting, the worst case is about 8s. Table ^ has a list of the periods used in 
this diagram (minus multiplets and repeated modes) along with the k assignment in the only i=l 
model. 
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The observed AP vs. P diagram for G29-38. 
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Table 3. Observed periods used in Figure ^ 
The k assignments assume all the modes are 
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Current published models (Bradley 1993, 1995) do not do a very good job of fitting this 
diagram, but we are encouraged by this as it means there is something we can learn from new 
modelling efforts. What is different about this diagram is the trend, or perhaps sudden change, 
to lower period spacing with higher period. This trend could be revealing physics unaccounted 
for in the model star, or it could mean some of these later modes are not i=l, as some of the 
alternative identifications plotted in Figure || suggest. Taking the dotted line path, this trend 
effectively disappears, but we have to come up with new assignments for the skipped-over modes. 
Bradley & Kleinman (1997) have shown, however, they can get a reasonable DA model to match 
the observed modes assuming mostly i=l. 

Based on the i=l only interpretation of the observed modes, we calculate a mean period 
spacing, P^, of 47 it 12s. This value can increase by a few (3-4) seconds, with the standard 
deviation decreasing an almost equal amount, if we take some of the alternative identifications. 
For now, we will work with the 47s spacing, determined by assuming all the modes are the i=l 
and have the k assignments made earlier. (We have here called the two modes near A;=17, k=17 
and k=18.) 

Using models of Bradley (1996) with a hydrogen layer mass of 1.0 x lO^^M^, consistent with 
the results of Clemens (1993, 1994), and a temperature of 11820 K from Bergeron et al. (1995), 
this period spacing corresponds to a mass of roughly O.GOMq, quite in line with the mean observed 
white dwarf mass. For this comparison, we used the 11700 K, O.6M0, standard model with 
1.0 X 10~^M^ hydrogen layer and helium layer of 1.0 x 10~^M-^., of Bradley (1996) with its mean 
i=l period spacing (determined from k=l to k=15) of 47 it 14 s. Unfortunately, this particular 
model does not match the individual modes very well, but it helps to establish the observed period 
spacing is in the correct range for a reasonable choice of parameters and to note the standard 
deviations of the model spacings and our observed spacing are similar. 

The mean period spacing is mainly a function of overall stellar mass. The hydrogen layer 
mass, however, also has an effect, and once we allow it to be a free parameter, we must consider its 
effects when determining a model asteroseismological fit. If we have another mass measurement, 
we can fix the overall stellar mass and just vary the hydrogen layer mass to match the observed 
mean period spacing. The most recent G29-38 mass estimate is from Bergeron et al. (1995). They 
derive a value near 0.69 Mq. We can fit this mass with our Pq if we use a 10~^''M^ hydrogen layer 
mass (iBradley 1993| ), although once again, the detailed mode list doesn't match very well. 

The most critical mismatch of the above-mentioned models to our list of modes are the 
k=l and k=2 modes. The O.6M0 model, for example, has its k=l, 1=1 mode at 145 s; the 
O.7M0 model has its k=l, i=l mode at 210s. Our labelled k=l and k=2 modes are at 110s 
and 177s respectively. There is clearly a mismatch here. The observed k=l and k=2 modes are 
certainly the lowest amplitude of the other modes, and therefore more prone to uncertainty, and 
could possibly be spurious. The case against this, however is fairly strong, as all the frequencies 
identified with lower amplitudes than these turned out to be linear combination frequencies; hence 



Tabic 3 ContiiiTiod 



Period 



18 915 
?? 1147 
?? 1240 



^This mode has a sum and difference 
combination with the k=10 mode, thus the 
"real" non-combination mode is either 110s 

as shown here, or 134s. 

^Thc identity of this mode is questionable. 
There is excess power in this region, but 
isolating it to a particular frequency is 
difficult. It should not be a strong constraint 
in model-fitting attempts. 



- 26 - 



their amplitudes are indeed significant. It could also be these modes are part of unidentified 
combinations, but there is as yet no evidence for this. It would be nice to see these modes appear 
again in additional (perhaps archival) data to help answer these questions. Given the ambiguities 
in the identification of these two modes already discussed, modelling efforts will have to be very 
careful and rigorous with all the possibilities here. 

There is one thing we can say with a much greater confidence: the identified set of 
non-combination peaks are not predominantly i=2 modes. Such a spacing with i=2 modes implies 
a stellar mass near 0.2M^, well out of the range of all previous mass estimates. The majority of 
the modes must be i=l. 

Preliminary results ( Kleinman 1995 , Dolez &: Kleinman 1997| ) show that where modes have 
been identified in other cool DAVs (from a slightly less comprehensive data set than that presented 
here), there are usually modes in G29-38 as well. This result extends that of Clemens to include 
the cooler DAVs and will be published in a future paper. If we follow the same line of reasoning 
as Clemens, this uniformity also suggests we have identified predominantly i=l modes. 



8. Conclusions 

After many years of searching, we finally have a DAV, G29-38, with enough observed modes 
to make a detailed asteroseismological analysis possible. With the extensive, 10-year data set 
presented here, we have separated the linear combination frequencies from the normal modes 
that recur in the star's Fourier spectrum. We have uncovered a discrete set of modes from 110s 
to 1193s with only a few gaps. We present evidence to support the assertion that these modes 
are predominantly i=l pulsations. There have never been so many modes in a DAV offered for 
analysis; we therefore expect surprises. Preliminary modelling efforts do indicate we can match 
these modes with current models, but so far the matches are not quite as close as we would prefer. 
These efforts are only preliminary as we have not yet attempted a systematic search of parameter 
space to fit these observations. 

Depending on the exact value of the Po, the period spacing, and the mass of the hydrogen 
layer, there may be a discrepancy between the asteroseismological mass and other measurements. 
The solution to this potential disagreement is as yet unknown, but perhaps more careful modelling 
efforts will guide us to the answer: are our mode identifications incorrect, the models lacking, 
or is G29~38's behavior affecting the other mass-determination methods, biasing them to an 
incorrect answer? Undiscussed since the Introduction, as it seemed to be completely independent 
to the presentation up to this point, however, is G29-38's infrared excess. Clearly something 
strange is happening to this star and until we definitively tie down the cause of the infrared 
excess, we cannot rule out its possible effects on the pulsations. We are encouraged, however, that 
Zuckerman (private communication) reports G29-38's IR excess remains unique among the DAVs 
and that so far, G29-38's pulsation spectrum looks similar to other cool DAVs. An intriguing 
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possibility, however, involves the effect of the large amplitude pulsations on an H-layer close to the 
nuclear-burning limit. Could this provide the source for both the pulsation instabilities and the IR 
excess? Koester et al. (1997) suggest not. They report atmospheric detection of heavy elements on 
the surface of G29-38, consistent with an increased accretion rate from surrounding circumstellar 
material, perhaps settling the mystery of G29-38's infrared excess. 

To help address these questions, we have obtained data on many of the other cool DA 
pulsators. Preliminary results do not yet address the mass problem (as not enough modes have 
been identified), but the modes we have seen all appear nearly identical to those in G29-38. This 
result, when added to Clemens's earlier work on the hotter DAVs, really says the DAs are a very 
homogeneous class. Therefore, given this set of modes observed in G29-38, it is extremely unlikely 
that they are anything but i=l pulsations such as have been found in the hotter DAVs. To 
summarize in broad terms, then, the evidence presented here for a set of i=l modes with a near 
47s spacing in G29-38 is the following: 



• There is a fairly regular ~47s period spacing observed. This value fits well in the range 
expected for stars with masses between 0.5 and O.7M0. 

• The mean period spacing is too high to be i=2 and maintain a reasonable mass for G29-38. 

• When we see multiplets, we see only triplets, consistent with only i=l rotationally-split 
modes. 

• The cool DAVs as a class share many common modes which appear to overlap those of the 
hotter DAVs which have already been suggested are £=1 modes ( [Kleinman" T995| and polej 



fc Kleinman 19971) . 



If we are going to label these modes as something other than i=l g-modes, we have to explain 
these four observations as coincidences. One accidental agreement is relatively easy to dismiss; 
four are much more difficult. There is no easy way to explain the majority of these modes as 
anything but i=l. 

The suggestion of global uniformity in the DA stars is so profound, we must look at it more 
closely. In particular, we must direct our efforts at determining both the stellar and hydrogen 
layer masses. We are now continuing to observe other DAs to determine the extent of their 
homogeneity. In February, 1996, the WET targeted the cool DAV, HL Tau 76. We hope data 
from this run, plus previous single-site runs, will uncover a similarly rich pulsation spectrum as 
we found here. Hopefully, the results form these observations will help guide us to some of the 
answers we posed earlier. These results should also help address the source of our model fitting 
problems, should they remain after a more detailed attempt. 

No matter what the numerical results of future modelling are, we have now made the 
necessary first step for their analysis: we have shown the cool DAVs are normal- mode pulsators 
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whose stable modes reappear in precise places in the power spectra, amidst a veritable forest 
of combination modes. We have taken a star whose power spectra seemed unsolvable, rapidly 
varying, and incapable of providing the kinds of clues needed for an asteroseismological analysis, 
and extracted such very clues. We must now, not only carefully explore model space to match 
these observations, but obtain similar data on other stars similar to G29-38 so we can begin to 
determine which conflicts between star and model belong only to G29-38, and which belong to 
the entire group of stars and therefore, to our understanding of the general physics which we put 
into the models. An increased understanding of the physical nature and evolution of the white 
dwarfs and the chance to calibrate more accurately the white dwarf cooling sequence, and hence 
the age of the local galactic disk, await only models to fit to these and future observations. 
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